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Abstract 

We present analytical and numerical calculations for some exited states 
of the frustrated J\—J% spin-| Heisenberg model for linear chains and 
square lattices. We consider the lowest eigenstates in the subspaces de- 
termined by the eigenvalue M of the spin operator Sl otal . Because of 
the reduced number of Ising basis states in the subspaces with higher M 
we are able to diagonalize systems with up to N = 144 spins. We find 
evidence that the Marshall-Peierls sign rule survives for a relatively large 
frustration parameter J2. 

PACS numbers: 75.10.J, 75.40.M 



1 Introduction 

We study the Marshall-Peierls sign rule (MPSR) for the frustrated spin-i Heisen- 
berg antiferromagnet 

H = Ji SiSj + J 2 £ SiSj, (1) 

(NN) (NNN) 
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where (NN) and (NNN) denote nearest-neighbor and next-nearest-neighbor 
bonds on a linear chain or a square lattice. 

According to Marshall's early work in 1955 we know the relative phases 
of the Ising basis states building the ground state wave function of a bipartite 
spin-i Heisenberg antiferromagnet (MPSR, for more details see below). Later 
on Lieb and Mattis generalized the theorem to arbitrary site spins and bipartite 
lattices without translational symmetries ^ ^ . The knowledge of the sign is of 
great importance for the construction of variational wave functions (see e.g. |Q) 
and for quantum Monte Carlo procedures which may have the so-called minus- 
sign problem if the MPSR is violated (see e.g.||). In particular, the possible 
violation of the MPSR in frustrated systems is a serious problem for variational 
and quantum Monte Carlo procedures. 

For the considered model (|l|) the MPSR can be proved for J2 < 0. In the 
frustrated model (J2 > 0) the MPSR can be destroyed {§ |, [|. However, in 
a recent work (J, Q we have presented general arguments that the MPSR may 
survive for relatively large Ji . These arguments are based on exact diagonaliza- 
tion results for small clusters (number of sites N < 24), as well as on the spin 
wave approximation. Following these arguments Zeng and Parkinson || use 
the MPSR as a new way of investigating the spatial periodicity in the ground 
state of frustrated spin chains. Furthermore, they studied the breakdown of the 
MPSR as a function of the chain length and the frustrating J2. By finite size 
extrapolation they estimated a finite critical value for Ji for an infinite chain of 
about 0.03Ji where the MPSR is violated in the ground state. 

Because of the exponential growth of the number of basis states the direct 
numerical calculation of the singlet ground state is limited to small clusters and 
the conclusions obtained from small systems seem not to be quite reliable. 

In this paper we exploit the observation that the MPSR holds not only for 
the singlet ground state but also for every lowest eigenstate in any subspace with 
higher quantum number M < y of the z-component of the total spin. In these 
subspaces the number of basis states is reduced and it is possible to diagonalize 
much larger systems. With this data an approximation to the thermodynamic 
limit is more reliable. Below we present data up to N < 144 and we can conclude 
that the MPSR survives indeed a finite frustration J2 at least for states with 
higher quantum numbers M. 

2 Marshall Peierls sign rule 

In the unfrustrated limit J2 = the lowest eigenstate of the Hamiltonian (|l|) in 
each subspace of fixed eigenvalue M of the spin operator Sf otal reads 

*m = £4 m V>, cL M) >o. (2) 
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Here the Ising states \m) are defined by 

|m) = (-l) SA_MA |mi) ® |m 2 ) ® • ■ • <S> |mjv) , (3) 

where |mi), i = 1, • • • , JV, are the eigenstates of the site spin operator Sf (— s$ < 
nii < s^, S A = Y,ieA s i' M A(B) = Y,teA(B) m ^ M = M A + Mb- The lattice 
consists of two equivalent sublattices A and B. Si = s, i = 1, • • • , N. are the site 
spins. The summations in Eq.(||) are restricted by the condition Yli=i m i = M- 
The wave function (||) is not only an eigenstate of the unfrustrated Hamiltonian 
(J2 = 0) and St otal but simultaneously of the square of the total spin S^ otal 
with quantum number S =| M \. Because > is valid for each m from 
the basis set (^|) it is impossible to build up other orthonormal states without 
using negative amplitudes Cm''. Hence the ground state wave function is 
nondegenerated. 

For the lowest energy eigenvalues E(S) belonging to the subspace M we 
have the Lieb-Mattis level-ordering 

E(S) <E(S + l) , S>0 , (4) 

i.e. the ground state is a singlet. Note that this level ordering might be violated 
by frustration. However, a lot of numerical calculations show the same level 
ordering also for strongly frustrated systems || |l(J . 



3 Analytic Solutions 

Now we include frustration (J 2 > 0). In the subspace with maximum M = N/2 
the MPSR is never violated in any dimension. Here the only possible state is 
the fully polarized ferromagnetic state which does not change with increasing 

In the next subspace M = (N/2) — 1 analytic solutions are found for linear 
chains and square lattices. In this subspace we deal with the so-called one- 
magnon state, where the wave function can be expressed as a Bloch-wave with 
a given k. 

Linear Chain - In this case the solution can be found by comparing the 
energies as a function of the wave vector k. 

N N 
E{k) = Ji(— - 1 + cos(fc)) + J 2 (— - 1 + cos(2fe)) (5) 

with k = = Q,±l,±2, . . .,+f. The comparison of E(w) and E(^(f -1)) 
yields the equation for the critical J 2 

i+ cos Hj-M (6) 

1 - cos [2n(l -I)] 
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In the limit N — > oo one obtains Jf = \ J\- 

Square lattice - For small J 2 in the considered subspace the lowest energy is 
obtained for k = (tt,w) and reads E\ = J%(N — 8) + J2N. The corresponding 
eigenfunction fulfills the MPSR. For larger J 2 we have to distinguish between 
two cases: (a) If the number of spins in the sublattice N/2 is even we find a 
transition at J 2 = ( Ji/2) to a twofold degenerated ground state with k — (tt, 0) 
or k = (0, 7r) with an energy E 2 — Ji (N — 4) + J 2 (N — 8). This state violates the 
MPSR, i.e. we have Jf = \ J\. Notice, that the eigenfunctions with k — (tt,0) 
or k — (0,7r) fulfill the so-called product-MPSR Qj. (b) If the number of spins 
in the sublattice is odd (e.g. N = 26), the situation is more complicated. The 
energy levels cross each other for Jf slightly greater than | J\ . 



4 Numerical Results 

In subspaces with lower quantum numbers M < (N/2) — 1 we cannot find simple 
expressions for the eigenvalues and eigenfunctions. Hence, we present exact 
diagonalization data for M < (N/2) — 2. Using a modified Lanczos procedure 
we calculate in every subspace M the state with the lowest energy E (M). Since 
the number of Ising basis states increases exponentially as ( N ^ M ), one can 
calculate E (M) for all M = ^ , . . . , only for relatively small systems (in our 
case N < 26). However, in subspaces with larger M we are able to present data 
for N up to 144. In all cases we use periodic boundary conditions. Analyzing 
the eigenfunction with respect to the MPSR we can determine numerically the 
critical J| where the MPSR is violated. 

Linear Chain - In Fig.|l] we show Jf as a function of l/N. For M(N) = 
(N/2) - 1 the analytic result is drawn. For the next M(N) = (N/2) - 2 the 
data show a similar behavior with the same critical value of J2 = \ for N — > 00. 
By considering the numerical data an analytic solution can be predicted 

J\ = ./, (7) 



1 + COS 


[ 7r ( 1 w=ry)J 


1 — COS 


2tt(1 - v £ r) )_ 



The following subspaces M(N) = (N/2) — p, p > 2 show a different behavior 
with different critical values for J 2 if N — > 00. These critical values decrease for 
increasing p but evidently a finite region with a non- violated MPSR is preserved. 

In Figj2| the critical Jf is shown as a function of a renormalized 
M r = M(N)/(N/2) for small systems (N = 8, . . . , 26) over the full range of 
M r . M r — 1 is the ground state subspace for a ferromagnet and M r — for an 
antiferromagnet. The finite size extrapolation for the ground state with M r = 
yields a small but finite critical value Jf w 0.03Ji which corresponds to the 
value estimated by Zeng and Parkinson in Q. The monotonic decrease of Jf 
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with decreasing M r indicates a finite region of a validity of the MPSR for all 

M r . 

Square lattice - In Fig.|| we show J| as a function of \/N. Here the N 
dependence is less regular and a finite size extrapolation is much more difficult. 
This is mainly connected with the shape of the periodic boundaries. For some 
of the finite lattices the boundaries are parallel to the x- and y-axis (e.g. for 
N =4x4, 6x6,. ..,12x12) whereas for other lattices the boundaries are inclined 
(e.g. N = 18, 20, 32, see e.g @). The impression of an oscillating behavior of J| 
versus 1 /N stems just from the alternation of parallel and inclined finite lattices. 
Nevertheless, it is evident that the critical Jf goes to a finite value for N — » oo. 
An extrapolation to the thermodynamic limit for the antiferromagnetic ground 
state, i.e. subspace M = 0, is almost impossible for the square lattice. However, 
if we assume for M = that the Jf is almost independent of 1/N for N > 16 
(as it is suggested by Fig.^J and by spin wave theory) we could estimate from 
our data for N = 10, 16, 18, 20 a critical value of about 0.2 . . . 0.3. 

Fig.[| supports this estimation. Here the critical J| is shown for different 
small lattices N < 34 as a function of M r . It is seen that for M r < 0.6 the 
critical Jf does not strongly depend on M r (in contrast to the linear chain, 
where we have a monotonic decrease) and gives for all the lattices a value of 
about 0.3Ji for the antiferromagnetic ground state (M r = 0). 

5 Conclusion 

For linear chains and square lattices we have shown, that in subspaces with 
large quantum number M of the spin operator S^ otal , the Marshall-Peierls sign 
rule is preserved up to a fairly large frustration parameter J£. 

Moreover, for linear chains the finite size extrapolation is quite reliable even 
for the singlet ground state and yields for all M a finite parameter region for 
J 2 where the MPSR is valid. 

For square lattices we observe in general higher critical values Jf compared 
to linear chains. From this observation and from the extrapolation for subspaces 
with higher M we argue that for square lattices the MPSR is stable against a 
finite frustration in all subspaces, too. 
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Figure 1: Linear chain: The critical value Jf where the MPSR breaks down 
versus 1/N for different M (J 1 = l). 



Figure 2: Linear chain: The critical value Jf where the MPSR breaks down 
versus M r for different N (Ji = 1). 



Figure 3: Square Lattice: The critical value Jf where the MPSR breaks down 
versus 1/N for different M (J 1 = l). 



Figure 4: Square Lattice: The critical value Jf where the MPSR breaks down 
versus M r for different N (Ji = f). 
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